Novel nonreciprocal acoustic effects in antiferromagnets 
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The possible occurrence of nonreciprocal acoustic effects in antiferromag- 
00 ' ne ts in the absence of an external magnetic field is investigated using both 

(i) a microscopic formulation of the magnetoelastic interaction between spins 
and phonons and (ii) symmetry arguments. We predict for certain antifer- 
romagnets the existence of two new nonreciprocal (non-time invariant) ef- 

o ' 

fects: A boundary-condition induced nonreciprocal effect and the occurrence 
of transversal phonon modes propagating in opposite directions having differ- 
ent velocities. Estimates are given and possible materials for these effects to 
be observed are suggested. 
PACS numbers: 72.55 +s, 73.50.Rb, 85.70.Ec 

Introduction It is well known in optics that the coupling of light to the order parameter 
in a ferromagnet mediates nonreciprocal effects such as the Faraday effect [|TJ . The Faraday 
effect in ferromagnets appears as a rotation of the plane of polarization of light which is 
incident along the axis of magnetization. An analogous effect in acoustics (the so-called 
acoustic Faraday effect) was first predicted by Kittel ||. Subsequently, it was verified 
experimentally |§ that in a ferromagnet, circularly polarized phonons of different sense have 
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different velocities when the direction of propagation is parallel to the axis of magnetization. 
Such an effect can be understood in terms of the magnetoelastic coupling between spins and 
phonons. Boiteux et al. Q performed in 1971 a detailed study of the magnetoacoustic 
effects in antiferromagnetic C^Os. They found that the crystal symmetry of Ci^Os allows 
an acoustic Faraday effect only in the presence of an external magnetic field. 

The nonreciprocal effects described above rely on the existence of a well defined magnetic 
moment induced either by the internal moments in the case of ferromagnets, or by the 
application of an external magnetic field in the case of C^Os ||. It is not apparent that 
such effects can occur in ordered antiferromagnets in the absence of an external magnetic field 
where the total magnetic moment is zero. In the case of optics, symmetry arguments can be 
used to deduce the presence of nonreciprocal bulk effects as well as a range of characteristic 
surface effects || in antiferromagnets. Indeed, nonreciprocal effects in antiferromagnets 
were demonstrated recently using the technique of Second Harmonic Generation 0. A 
microscopic explanation of these experiments which relies on the symmetry of the crystal 
and the evaluation of the interactions at the atomic level has also been proposed ||[| . 

In this Letter we examine if it is possible to observe nonreciprocal phonons in antiferro- 
magnets in the absence of an external magnetic field. We find that the answer is affirmative 
but these effects are very different from those observed in optics. Our calculations show that 
there is a new boundary-condition induced elliptical-polarization effect in antiferromagnets. 
In addition we predict, for certain antiferromagnets, the occurrence of transversal phonon 
modes propagating in opposite directions having different velocities, an effect which cannot 
occur in ferromagnets. 

Ferromagnets Kittel was the first to predict the magnetoacoustic Faraday effect in 
ferromagnets within a semiclassical continuum formulation. For the sake of completeness, 
we present a brief derivation of Kittel's result within a linear chain model for a ferromagnet. 
We take the z-axis to be the direction of the chain of ions corresponding to the main 
crystallographic axis, which we assume to have C4 or higher rotational symmetry. The 
magnetization is parallel to the £-axis. 
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We define u na and p nn (n denotes the site index and 7 = x, y) to be the canonical 
displacement and momentum operators obeying the commutation relations [w„ )7 , Pn'rf] = 
ih8 ntn i8y t yi- The Hamiltonian H = T + V e \ + Hneis + contains the unperturbed phonon 
part 

T = ^T,P 2 n,V V * = TEK+LT-WnJ (1) 

where K is the elastic constant, and m the mass of the ion. Hneis corresponds to the 
Heisenberg Hamiltonian and describes the unperturbed spin-spin interaction, 

Hueis = J*^2S n ■ S n+ i (2) 

n 

J < is the ferromagnetic exchange constant and S n corresponds to the spin operator 
on site n, with [S%, S^J = ih8 n , m e a ^S1 (a, /3, 7 = x, y, z). H sp describes the magnetoelastic 
coupling between phonons and spins: 

H sp —9^2 ( M n+1,7 — u n,j) ('S'n+l'S'n + ^n+l^n) ■ (3) 

g is the magnetoelastic coupling constant. 

Eq. (^) is rotational invariant around the i-axis. A 90-degree rotation around the z- 
axis would correspond to (u n ,x — > —u njy , u n , y — > u njX , — > —S%, S% — > Eq. @ is 

the lattice version of the usual coupling between components of the strain tensor and the 
magnetization M . This interaction term can also be understood as a higher order expansion 



of the exchange term and has its origin in the spin-orbit coupling fllO|| . In order to obtain 
the time dependence of the variables p nn , u nn and S% we consider the Heisenberg equation 
of motion ih-^A = [A, H\. For simplicity, we work with circular-polarized coordinates = 
u U)X ± iu nt y, S„ — S® ± iS% and we assume a plane wave ansatz u^(t) = exp[i(ut — nk)]^, 
S^(t) = exp[i(ut — nk^S"* 1 (the unit cell constant is taken to be one) in order to solve the 
equations of motion. We set the magnetization to be along the z-axis ((S*) = 0, (S^) = 0, 
(S z ) = (S z )) and we obtain after eliminating p nn : 

uS ± = ±g(S z ) 2 {2t sin tyu* ± J{S Z )2{\ - cos k)S ± (4) 
(-muj 2 + 2K(1 - cosA;))-u ± = g(S z )(—2isink)S ± . (5) 



Combining Eq. (|j) and (||) we get: 

{[-mu 2 + 2K(1 - cosA;)][w T J(S Z )2(1 - cos A;)] T V(5 2 ) 3 (sin A;) 2 } m ± = (6) 

which describes the dispersion relation of the phonon system interacting with the spin sys- 
tem. Note that for u + and u~ the dispersion relation is different implying that for a given 
frequency lu, right (+) and left (-) circular polarized phonons have different wave vectors k + 
and k~ respectively. 

An incoming linear polarized wave with frequency u is decomposed into a linear super- 
position of right and left circular polarized waves with the same frequency uo but different 
wave vectors k + /k~ via exp[i(ujt — k + n)] u + + exp[i(ut — k~n)\ u ~. The physical amplitudes, 
determined by the real part of u are given by 

u n ,x = ucos(ut — n{k + + k~)/2) cos(n(k~ — k + )/2) 



u 



n,y 



ucos(ut - n{k + + k~)/2) sm(n(k~ - k + )/2). 



The plane of polarization rotates around the z-axis and it is easy to obtain from the 
solution of Eq. @ the rotatory power {k + — k~)/{k + + k~) « g 2 (S z ) 3 /(Kuj). The nonre- 
ciprocal behavior is a direct consequence of the coupling between phonons and spins in a 
ferromagnet. For a typical ferrite the rotatory power should be of the order of 1/4 for sound 
waves in the microwave regime (w ~ 10 9 s~l), as was corroborated experimentally ||. 

Antiferromagnets As a model for an antiferromagnetic material, we consider a linear 
chain with two magnetic ions {A/B) per unit-cell, with UA, n and ub, u being their respective 
displacement vectors. A straightforward generalization of ([!]), (|2]) and (^j leads to 



2m „ )7 

V e l = XX M -B,n,7 — u A,n,j) 2 + XX U An+l,7 _ U B,n,~f) 2 • (7) 

n,7 ^ n,7 

for the phonon part. K\ and K2 are, respectively, the intra- and intercell elastic constants. 

Hueis = [^l^A" ' £>B,n + ^'S's.n " <?A,n+l] > (8) 



corresponds to the Heisenberg term with Ji > and J 2 > (antiferromagnetic). And the 
magnetoelastic coupling is given by 



Hsp — 9l y^X u B,n,-y u A,n,f)(S'^ n S An + ^B,n^A,n) 
+02 y^(^A,n+l,7 ~~ U B,n^){S\ n+ iS Z Bn + ^n+l^n) 



(9) 
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where gi and 52 are, respectively, intra- and intercell magnetoelastic coupling constants. In 
order to obtain the dispersion relation of the phonons interacting with the spin system, we set 
the staggered magnetization to be parallel to the z axis ({S An ) — > (S A ) and (Sg n ) — > — (S A )) 
and solve the equations of motion for ua/b,u, Pa/b,u and Sa/b,u- 

We consider the plane- wave ansatz u n = exp[i(ut — kn)} u and define the vector u n = 
(uA,n,x, ub,u,x] u A ,n, y , u B ,n, y )- Eliminating p A /B,n,-y and S A /B,n,"f we obtain Mu = 0, where the 
4x4 hermitian matrix M is given by: 



L(u) + 2Js A R k -T k (u) - Js 4 P k 2is 3 Q k uj 
-T_ fc (cj) - Js 4 P_ fe L(uj) + 2Js i R k 

-2is 3 Q k uj L(u) + 2Js 4 R k 





— 2is 3 Q k u 
-Tduo) - Js A P k 



\ 



V 



2is 3 Q k u -T_ k (uu) - Js 4 P_ fc L(u) + 2Js 4 R k 



s stands for (S A ), uj 2 = 4 J 2 s 2 sin 2 (k/2). We consider for simplicity J\ = J 2 = J and the 
following definitions are used: 

Qk = 0i02(l - cos/c) 
= (0i + 02 ~ %i0 2 )(l + cos k) + 2(# 2 + g 2 2 ) 
P k = (1 + exp[-ifc])(#i - 5f 2 exp[iA;]) 2 + 
(1 + exp[ifc])(^i - 2 ) 2 + 4(^1 - 5- 2 )(5 , i - 02 exp[?A;]) 
= (-mw 2 + Ki + K 2 )(uj 2 - uj 2 g ) 
T k {uj) = {K x + K 2 exp[ik])(u 2 - u 2 ) (10) 



The eigenvectors of M are given by: 
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u + = (u A ,u B ;iu A ,iu B ) (11) 
u = (u* B , u A ; -iu* B , -iu* A ) (12) 

and the amplitudes u A and u B are determined by 

{L(uj) - 2us i Q k + 2Js A R k }u A = 

{T k (u) + Js 4 P k }u B (13) 

The dispersion relation uj = uj(k) is given by the solution of detM = 0. u + and u~ have the 
same dispersion what translates to the fact that right and left circular polarized phonons are 
degenerate and one might think, that there are no nonreciprocal effects in antiferromagnets. 
This degeneracy is broken when an external magnetic field is applied (see @]). We shall 
show here that, in fact, even in the absence of an external magnetic field, nonreciprocal 
effects might be obtained by considering an appropriate choice of boundary conditions. 

Let us consider a system that is driven by a linear polarized amplitude in the ^-direction 
at site n = 0; i.e., we apply the external driving force 

u A ,a,x = exp[iut] u {ext) , u A ,o, y = 0. (14) 
The appropriate superposition of right and left circular polarized acoustic waves u + /u~ that 



satisfy ( |14D is 

u = u* B u + + u A u' = (2u* B u A , \u B \ 2 + \u A \ 2 ; 0, i(\u B \ 2 - \u A \ 2 )), 

with 2u* B u A = u^ ext \ The physical amplitudes, given by the real part of exp[i(tot — kn))u, 
are therefore linearly polarized for the A-site, 

UA,n,x = COs(ut - kn) U {eXt \ U A ,n,y = 0, 

and elliptically polarized for the -B-sites: 

UB,n,x = cos(ut - kn) (\u A \ 2 + \u B \ 2 ), 
UB,n, y = sm(ujt - kn) (\u A \ 2 - \u B \ 2 ). 



The degree of the homogeneous elliptical polarization {\ua\ 2 — \ub\ 2 ) / (\ua\ 2 + \ub\ 2 ) ~ 
(S^) 3 gig2(l — cosk) / (mw 3 ) is proportional to the spin-phonon couplings gi, g 2 and to (S A ), 
which makes it nonreciprocal. Therefore, by choosing appropriate boundary conditions it 
is possible to observe this nonreciprocal effect in an antiferromagnet. Note, in this case, 
the nonreciprocity relates to the amplitudes of the sound wave and not to the velocity as it 
happens in the ordinary acoustic Faraday effect in ferromagnets. 

Since the magnetoelastic coupling between phonons and spins is determined by general 
symmetry arguments (see also H and J3J) that are not constrained by the dimensionality, 
we can extend our analysis to two and three dimensional materials. Then, it follows that 
the boundary-condition induced nonreciprocal effect should occur in materials that consist 
of antiferromagnetically stacked chains or planes with each chain/plane having a non van- 
ishing moment along the stacking direction. Typical examples are compounds belonging to 



the family of ABX3-type hexagonal antiferromagnets, like CSC0CI3 and CsNiBr3 ||1 1|| . The 
magnitude of the effect should be comparable to the nonreciprocal effect observed in ferro- 
magnets since both effects originate from the same microscopic coupling mechanism. 

Novel effects In the case illustrated above, we considered a system of ions aligned along 
the f-axis with cubic symmetry in the paramagnetic phase. Let us now assume a crystal 
with lower symmetry than cubic which also orders antiferromagnetically. We add to our 
Hamiltonian (||) the following term: 

H-sp 1 = fl y^X u B,n,x ~ u A,n,x){S V Bn S Z A n + S Bn S A ^ n ) 
n 

~fl y^X u B,n,v ~ u A,n,y){.S Bn S Z A n + S Bn S A ^ n ) 
n 

+ f2 XX M ArH-M — U B,n,x){.S V A n+1 S B n + S A n+l S B ^ n ) 
n 

—f% /XuA,n+l,y ~ u B,n,y){S X A)n+ iS Z B n + S An+1 S B ^ n ) (15) 
n 

where fi and f 2 are, respectively, intra- and intercell magnetoelastic coupling constants. 
Note that this term is invariant under the following symmetries: (1, 1, 2 Z , 2 Z , ±4 Z , ±4 Z ) but 
not, for instance, under a reflection 2_ xy . Such a term (H^') can be written for systems 
crystallizing, in the paramagnetic phase, in one of the following point groups: C 2 h, C^h, C31 



or C§h (in Schonflies notation). 

We calculate the dispersion relation for H = T + V e i + H Heis + + HjJ/ by considering 
the equations of motion for ua/b,u, PA/B,n and Sa/b,u- The staggered magnetization axis 
is taken to be parallel to the z-axis ((S A ^ n ) — ► (S A ) and (S% n ) — > — (S A )). The following 
characteristic polynomial is obtained 

[mV - 2muj 2 {K 1 + K 2 ) + 2K X K 2 {\ - cosA;)][c«; 2 - uj 2 } 2 

-4m(^) 3 a[sin k}{u 2 - oj 2 }uj 3 
-A{S z A f[gl + f 2 ][g 2 2 + f 2 ][l - cos*]V 
-4mJ(S z A } 4 [(B(3 + cos A;) - 4/3(1 + cosfc)]^ 2 - uj 2 ]^ 2 
+4(u 2 - u 2 ) J(S Z A ) 4 [3 + cos k] [1 - cos k] [(gf + f 2 )K 2 + (gj + f 2 )^} [u 2 - uj 2 } 

+ 8J 2 (S A n 9 2 1 + fM + f 2 }[l - cosA;] 3 = (16) 

where a = g 2 fi - f 2 gi, (3 = g x g 2 + fxf 2 , B = g\ + f 2 + g\ + f 2 , u is the dispersion relation 
of the unperturbed phonon system and J\ = J 2 = J. While the same dispersion relation 
uj = uj{k) is obtained whether we consider right (+) or left (-) circular polarized phonons, 
Eq. ([H]) has some important features. It breaks time reversal symmetry (0) u(k) ^ u(—k) 
due to the terms proportional to a sin k, and also inversion symmetry (I) but respects the 
product JO. The breaking of these two symmetries gives rise to a novel nonreciprocal 
effect entirely different from the Faraday effect. It translates to the fact that transversal 
phonon modes propagating in opposite directions have different velocities. Note that this 
nonreciprocal effect is directly proportional to a. If f\ and f 2 are zero, there is no such 
effect. 

^From a closer inspection of Eq. (|15D, it follows that the three conditions necessary for 
the above effect to occur in an antiferromagnet are (i) the point group should be one of 
the following: C 2 h C^, C 3i or C Qh fl^fl , (ii) the magnetic ions themselves are not centers 
of inversion and (iii) the magnetic moment has a component parallel to the main crystallo- 
graphy axis, (S z ) ^ 0. The wolframites P| CoW0 4 , NiW0 4 , FeW0 4 and FeNb0 4 fulfill all 
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three conditions and are therefore suitable candidates. Their crystal group is C2/1 and the 
magnetic unit cell (2a,b,c) is doubled with respect to the chemical unit cell. Another group 



of candidates are the ilmenites |T^] FeTiOs, NiTiOs and MnTiOs that belong to the point 
group C 3i . 

We estimate for MnTiOs, the relative difference 

u(k)-u(-k) 

u(k)+u(-k) { ' 

for the acoustic phonons. We obtain u(k) and u(—k) by solving Eq. fll6|) numerically. The 
values of the constants appearing in Eq. ([16]) are taken from the literature (see and 
|16| ). The magnetoelastic coupling constants g%, g-z, fi, J% for manganese ions in MnTiOs 
are not available in the literature. They can be estimated from first principles by doing an 



expansion of the exchange term up to the required order |10| . Here, for g 1 we assume values 
similar to those corresponding to other transition metals like Fe. We assume fi and fz to 
be somewhat smaller than Q\. This is reasonable given that they define a lower symmetry 
coupling than g\. MnTiOs has a spin- wave gap A(H) which varies strongly in an external 



magnetic field H ||15|| . We find the novel nonreciprocal effect to be maximal at the resonance 
frequency of the magnon and the acoustic phonon system which occurs at ui ~ 156 GHz. At 
resonance the ratio Eq. fll7|) is ~ 10 -4 . Note that this nonreciprocal effect occurs already in 
absence of an external magnetic field. In an external magnetic field the gap A (if) decreases 
linearly and vanishes at the spin-flop transition H c = 5.8 T. Since the resonance frequency 
is proportional to the spin-wave gap, a promising experimental set-up in order to measure 
Eq. (|TTD would be to sweep the external magnetic field for a fixed phonon frequency. Note 
that standard measurements of the velocity of acoustic phonons are accurate to 10 -5 — 10~ 6 . 

Conclusions We have studied the possible occurrence of an acoustic Faraday effect in 
antiferromagnets in the absence of an external magnetic field by using a Hamiltonian for- 
mulation of the magnetoelastic coupling between phonons and spins and general symmetry 
arguments. We predict a new boundary- condition induced homogeneous elliptical polar- 
ization effect for layered materials like the ABX 3 -type hexagonal antiferromagnets. The 



microscopic analysis of this effect shows that it should be of the same order of magnitude 
as the acoustic Faraday effect already observed in ferromagnets 0. We also predict the 
occurrence of transversal phonon modes propagating in opposite directions having different 
velocities in certain low-symmetry magnets like the wolframites or the ilmenites. An esti- 
mate of the order of magnitude of this effect for MnTi03 is given. We suggest experiments 
be undertaken to detect these effects. 
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